Nonlinear response of single-molecule nanomagnets: equilibrium and dynamical 
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We present an experimental study of the nonlinear susceptibility of Mni2 single-molecule magnets. 
We investigate both their thermal-equilibrium and dynamical nonlinear responses. The equilibrium 
results show the sensitivity of the nonlinear susceptibility to the magnetic anisotropy, which is nearly 
absent in the linear response for axes distributed at random. The nonlinear dynamic response of 
Mni2 was recently found to be very large and displaying peaks reversed with respect to classical 
superparamagnets [F. Luis et ai, Phys. Rev. Lett. 92, 107201 (2004)]. Here we corroborate the 
proposed explanation — strong field dependence of the relaxation rate due to the detuning of tunnel 
energy levels. This is done by studying the orientational dependence of the nonlinear susceptibility, 
which permits to isolate the quantum detuning contribution. Besides, from the analysis of the 
longitudinal and transverse contributions we estimate a bound for the decoherence time due to the 
coupling to the phonon bath. 

PACS numbers: 75.50.Xx, 75.50.Tt, 75.45.+j, 75.40.Gb 



I. INTRODUCTION 

Superparamagnets are nanoscale solids or clusters with 
a large net spin (S ~ 10 1 -10 4 ). This spin is coupled to 
the environmental degrees of freedom of the host mate- 
rial, e.g., phonons, nuclear spins, or conducting electrons. 
Due to the dynamical disturbances of the surroundings 
the spin may, among other things, undergo a Brownian- 
type "reversal" , overcoming the potential barriers cre- 
ated by the magnetic anisotropy. 

Depending on the relation between the reversal time 
r and the observation time t bs, different phenomenolo- 
gies can be found. For r <C i b s , the spin exhibits the 
thermal-equilibrium distribution of orientations as in a 
paramagnet; the large values of S are the reason for the 
name superparamagnetism. When r 3> t ^ s , in contrast, 
the reversal mechanisms appear blocked and the spin 
stays close to an energy minimum (stable magnetization 
conditions appropriate for magnetic storage). Finally, 
under intermediate conditions (r ~ i bs) one finds non- 
equilibrium phenomena (i.e., magnetic "relaxation"). 

For large S, for instance in magnetic nanoparticles , 
a classical description is adequate . The essential phys- 
ical ingredients are the thermo-activation over the mag- 
netic anisotropy barriers and the (damped) spin preces- 
sion @. As the spin value is reduced, quantum effects 
can start to play a role. For moderate spins (S ~ 10), 
as in single-molecule magnets Q , the quantum nature of 
the system comes significantly to the fore. For instance, 
the spin reversal may also occur by tunneling whenever 
the magnetic field brings into resonance quantum states 
located at the sides of the barrier (Fig. . 

Several fundamental problems can be studied on these 
systems 0,E1I3- First, the quantum-to-classical transi- 
tion, with the emergence of classical properties. Single- 
molecule magnets constitute a model system to study 
quantum mechanics at a mesoscopic level, while mag- 
netic nanoparticles provide a natural classical limit. Sec- 



ond, one can address the effects of environmental degrees 
of freedom on a given system. Classically, one faces the 
rich phenomenology of rotational Brownian motion of the 
nanoparticle magnetization. In the quantum case the 
bath coupling not only produces fluctuations and dissi- 
pation (allowing the system to relax to thermal equilib- 
rium), but it is also responsible for the decoherence of its 
quantum dynamics. Thus, single-molecule magnets con- 
stitute an important experimental benchmark to test the 
predictions of the theory of quantum dissipative systems 
[3 (much as Josephson junctions in superconductivity). 

The best studied magnetic molecular clusters are those 
named Mni 2 and Fes, both with S = 10 in the ground 
multiplet (for other examples see Ref. Q)- To under- 
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FIG. 1: Energy levels for Mni2 (the molecule is sketched in 
the inset along with the spin orientations of the Mn ions). 
The levels are plotted vs. the quantum number m at H = 
and show the bistable potential for the spin due to the 
magnetic anisotropy. The horizontal line marks the border 
between "classical" or localized energy levels, A m < £ m , and 
the tunneling levels A m > £ m (A is the tunnel splitting and 
£ the width of the environmental bias-field distribution) . 
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stand their behavior a plethora of theoretical calculations 
and every conceivable experimental technique have been 
used. However, the nonlinear susceptibi lity Y3, fruitfully 
exploited in studies of spin glasses 0, 0, [ljj, random 
anisotro py syste ms [12j. but also in magnetic nanoparti- 
cles [l3l \l4. Ilfij. had been overlooked. For classical su- 
perparamagnets X3 provides information on parameters 
to which the linear susceptibility is less sensitive, like the 
anisotropy constant D [l6j or the spin-bath coupling pa- 
rameter A (which enters the scene due to the strong 
dependence of the relaxation rate r = 1/r on transverse 
magnetic fields 18]). Besides, the dynamical nonlinear 
susceptibility has a genuinely quantum contribution due 
to the detuning of the energy levels by a longitudinal 
field 0]. It has a sign opposite to the classical (preces- 
sional) contribution, thus allowing to ascertain whether 
quantum effects, such as resonant tunneling, are relevant 
in a given nanomagnet (an issue sometimes controversial 

[13). 

In this article we present experimental results for the 
thermal-equilibrium and dynamical nonlinear suscepti- 
bility of Mni2 acetate. Compared to the linear suscep- 
tibility, the equilibrium ^3 shows an enhanced sensitiv- 
ity to the magnetic anisotropy, even for axes distributed 
at random (allowing to estimate D from measurements 
in powdered samples). As for the frequency-dependent 
X3, we study its dependence on the angle of the ap- 
plied field. This gives direct access to the relaxation-rate 
field-expansion coefficients, T — T\h=o oc g\\H^ + g±H]_, 
which contain valuable information on the spin rever- 
sal mechanisms, thus allowing to separate the "classical- 
transverse" and "quantum-longitudinal" contributions to 

X3- 

In the discussion simple approximate formulas and nu- 
merical results from the solution of a Pauli quantum- 
master equation are used. Our investigation confirms the 
interpretation of Ref. fl9| of the large quantum contribu- 
tion to X3(w) as arising from the detuning of the tun- 
nel channels by a longitudinal magnetic field. Thus, the 
experimental nonlinear resp onse is consistent with the 
established scenario |2l|, |22j, |23J, |2J, |2jJ |26( of thermally- 
activated tunnel via excited states in Mni2. The analysis 
also gives a bound for the decoherence time r$ (timescale 
for the attainment of a diagonal density matrix due to the 
coupling to the phonon bath). The obtained r$ turns out 
to be much shorter than the lifetime of the spin levels To, 
and is responsible for a fast loss of coherent dynamics 
(like tunnel oscillations or precession). 

II. SAMPLES AND MEASUREMENTS 
A. Samples and set up 

Single crystals of Mni2 acetate were grown following a 
procedure similar to that described by Lis [2I} . The con- 
centrations of the reactants, however, were higher than 
those of in order to increase the supersaturation and 
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FIG. 2: Magnetic ac susceptibility of a single crystal of Mni2, 
normalised by its zero-field value, vs. the static field H (par- 
allel to the anisotropy axes, ip = 0). Results for the real part 
at T — 5K and various frequencies lot are shown. The solid 
lines represent polynomial fits from which \3 is obtained. The 
parabolic approximation xi + 3%3H 2 , which dominates the 
low-field behavior (\H\ < 100 Oe), is shown by dotted lines. 

the growth rate, yielding larger crystals. These were re- 
grown several times by renewing the mother solution. X- 
ray diffraction patterns of powdered crystals agreed with 
simulated patterns from the known crystal structure. 

The measurements were performed on a single crystal 
with dimensions 3x0.5x0.5 mm 3 at different orientations 
with respect to the applied magnetic field. To this end, 
we constructed a rotating sample holder that enables the 
c crystallographic axis (which defines the anisotropy axes 
of the Mni2 molecules) to be rotated a given angle tp with 
respect to the magnet axis. This angle is measured with 
a precision better than 0.5 degrees. To calibrate the posi- 
tion of the zero we used the measured linear equilibrium 
susceptibility, which should be maximum when the field 
is parallel to the anisotropy axis (ip = Q). 



B. Magnetic measurements 

Dynamical susceptibility measurements were per- 
formed using the ac option of a commercial SQUID mag- 
netometer, by applying an alternating field ~ Ahe 1U3t . 
The ac susceptibility was measured under a weak super- 
imposed dc field H, parallel to the oscillating one. The 
first harmonic of the response can then be expanded as 

H) = xiM + 3 X3 M# 2 + 5 X5 M# 4 + • ■ ■ • (i) 

The (^-independent) expansion coefficients give the or- 
dinary linear susceptibility xi an d the nonlinear ones X3j 
X5,- ■ ■ ■ As it is customary, we focus on X3 and refer to it 
as the nonlinear susceptibility. 
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FIG. 3: Magnetic susceptibility vs. frequency measured along 
the anisotropy axis (ip = 0) at T = 5K. Results are shown at 
zero bias field (circles) and at H = 300 Oe (squares). Solid 
symbols are for the real part (in-phase component) and open 
symbols for the imaginary part (out-of-phase) . The lines are 
fits to the Debye law (|5J with t\h=o = 1.3(3) x 10~ 2 s and 
t|h=300 = 2.8(1) x HT 2 s. 

To determine the nonlinear susceptibility we per- 
formed polynomial fits of the ^/-dependent ac data whose 
quadratic coefficient gives \3- An illustrative example 
of the fitting procedure is shown in Fig. [5J For suffi- 
ciently low H, a good description is provided by a simple 
parabolic dependence Xi+3x3-ff 2 - For increasingly larger 
fields, we increased the order of the polynomials when- 
ever the fitting error became greater than 5%. The exper- 
imental X3 was taken as the mean value of all quadratic 
coefficients obtained from the different order polynomi- 
als, thus miminizing the error of the determination. 

The measurements were performed at temperatures in 
the range 2K < T < 45 K and frequencies 0.08 Hz < 
ui/2tt < 1.5 kHz. The amplitude of the ac field (Ah = 
3 Oe) was sufficiently small not to induce nonlinearity 
in Ah (associated to the generation of harmonics). Ac- 
tually, the X3( w ) obtained here from the field-dependent 
susceptibility differs from the X3 extracted from the third 
harmonic of the response to an ac field (at H = 0) . How- 
ever, as we show in Appendix^ the dependence of both 
quantities on the parameters under study is analogue. 



C. Equilibrium vs. dynamical measurements 

In the next sections we are going to study the equilib- 
rium and dynamical susceptibilities. Let us define here 
a practical criterion to decide when the experimentally 
measured Xi(w) and X3( w ) correspond to equilibrium or 
off-equilibrium conditions. In the temperature range cov- 
ered by our experiments, T > 2K, the molecular spins 
of Mni2 relax via a thermally-activated tunneling mecha- 
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FIG. 4: Linear (top panel) and nonlinear susceptibilities (bot- 
tom) vs. temperature for several frequencies (only 9 Hz for 
X3) measured along the anisotropy axis. Solid symbols are 
for the real parts and open symbols for the imaginary parts. 
Dashed lines are the equilibrium susceptibilities in the Ising 
(large anisotropy) limit [Eqs. © and Hul l. The vertical line 
marks the boundary above which thermal equilibrium results 
are safely obtained using u>/2tt — 9 Hz. Inset: temperature 
dependence of the reciprocal equilibrium linear susceptibility 
and its fit to a Curie- Weiss law |g) for T < 10 K. 

nism 0,113, [2^. This process gives rise to a well-defined 
relaxation time r and the ac response can be described 
by a simple Debye formula: 



X = XS + 



XT - Xs 

1 



1UT 



(2) 



Here xt and xs are the isothermal (thermal equilibrium) 
and adiabatic limits of x- In Fig. [3] we show how the 
response of the Mni2 crystals follows Eq. (|2J at tempera- 
tures and magnetic fields typical of our experiments |28j . 
The equilibrium regime corresponds to frequencies fulfill- 
ing ujt <C 1 (left part of the plot), relaxation effects be- 
coming important when the range lot ~ 1 is approached. 

On the other hand, the relaxation time of this system 
increases exponentially as T decreases, following an Ar- 
rhenius law (see for instance Fig. 2 in Ref. [Hj) 



t = T exp (U/ksT) 



(3) 
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Here U is an activation energy and To an attempt time, 
which set the magnitude and temperature dependence 
of r. At zero field we obtained Uq ~ 65 K and r ~ 
3 x lCT 8 s 19]. In a x vs. T experiment (Fig. the 
condition lot — 1 defines a superparamagnetic "blocking" 
temperature k B Th = — Uq ln(wTo). Below Tb, the real 
part x'i drops from \t towards xs whereas x'i departs 
from zero and shows a maximum. As for the nonlinear 
susceptibility, the non-equilibrium response near Tb leads 
also to a non-zero imaginary part x'ii and to a strong 
deviation of x'3 from the equilibrium X3T- in contrast 
with the linear response (which decreases from Xt)i x'3 
becomes larger than X3T near Tb 0|i this dynamical 
nonlinear phenomenon will be discussed in detail in the 
following sections. 

As can be seen in Figs. and 0J the transition from 
isothermal to adiabatic conditions extends over a cer- 
tain frequency or temperature range, determined by the 
width of the x"( w i T) curve. As a practical rule, we take 
Xt = Xi an d X3T — x'3 when the imaginary parts are rea- 
sonably small x"/x' < 10~ 2 . This value is, as a matter of 
fact, close to the experimental error of the present mea- 
surements. Using this criterion we have extracted, from 
the dynamical susceptibility data vs. T, the equilibrium 
Xt and X3T discussed in the next section. 



III. EQUILIBRIUM LINEAR AND 
NONLINEAR SUSCEPTIBILITIES 



The equilibrium xt obtained as described in the pre- 
vious section is shown in the inset of Fig. 21 In the tem- 
perature range 4K < T < 10 K it approximately follows 
a Curie- Weiss law 



Brillouin magnetization ( Appendix I A lfl . and read 



N A 



X3iso = -Njl 



(g f i B ) 2 S(S+l) 
3k B T 

45(fc B T) 3 



(5) 



s(s + i)[s(s + i) + h] , (6) 



where is the number of molecules per mol. The first 
of these is merely Curie's law and the second its nonlin- 
ear counterpart. Normalized by Xiso and X3iso, the ex- 
perimental susceptibilities lose their bare 1/T and 1/T 3 
contributions, so that their remaining temperature de- 
pendence is mostly due to the effects of the anisotropy. 

The so normalized equilibrium susceptibility data are 
shown in Fig. [S] Clearly, the isotropic limit is only at- 
tained for sufficiently high temperatures (T > 30 K). 
Note that the high-T limits of xt and X3T are slightly 
smaller than Xiso and X3iso- This is caused by the ther- 
mal population of higher-energy spin multiplets of the 
cluster, the lowest of which has S = 9 in Mni2. There- 
fore, in that temperature range the Mni2 molecule can no 
longer be seen as a superparamagnetic spin S = 10 and 
the thermal mixture of spin states reduces its suscepti- 
bilities (the analogue to the decrease of the spontaneous 
magnetization in a solid by excitation of spin- waves). 

As the temperature decreases, both Xr/Xiso and 
X3r/X3iso increase, departing from ~ 1. This is natural 
since Eqs. iJSJ-JfjJ are only valid when the thermal en- 
ergy k B T is larger than all zero-field splittings (produced 
by the magnetic anisotropy). The simplest Hamiltonian 
that describes the magnetic behavior of an isolated Mni2 
molecule contains the Zeeman plus uniaxial anisotropy 
terms (see also Fig. ^ 
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with a Curie temperature 9 ~ 1.2(2) K. Measurements 
performed on a powdered sample (not shown) , which are 
less affected by anisotropy effects, also follow Eq. Q al- 
beit with a smaller 9 = 0.5(1) K [for the interplay of inter- 
actions and anisotropy, see Eqs. (3.2)-(3.6) in Ref. p9j]- 
These finite 9 point actually to the presence of dipolar 
interactions between the molecular spins in the crystal, 
which would give rise to long-range order at sufficiently 
low temperatures. In the (super)paramagnetic regime of 
interest here, interactions merely produce a susceptibility 
somewhat larger than that of non-interacting clusters. 

In addition to interactions, a most important influence 
on the temperature-dependent susceptibilities is exerted 
by the magnetic anisotropy [3(]]. To illustrate these ef- 
fects, it is helpful to normalize the experimental xt and 
X3T by their isotropic limits, Xiso and X3iso- These can 
be obtained from the field expansion coefficients of the 



Here D ~ 0.6 K and A4 ~ 10 _3 K are the second and 
fourth-order anisotropy constants for Mn 12 |31|. and 
H x ,y,z the components of the field along the (a, 6, c) 
crystallographic axes. The largest zero-field splitting 
produced by the anisotropy occurs between the states 
to = ±(5 — 1) and the ground state to = ±5: 

ft - (2S 1 - 1)D + [S 4 - (S - \f]Ai ~ 14.8 K . (8) 

When k B T becomes comparable to Slo several related ef- 
fects occur: (i) the magnetization is no longer given by 
the Brillouin law and T does not appear in the combi- 
nation H/T, (ii) xt and X3T deviate from the simple 
Eqs. ©-© and depend on tp, and (iii) the normalised 
susceptibilities adquire a dependence on T. For classical 
spins, these effects were studied in Ref. |32l |. 

Eventually, in the low temperature limit k B T/D — > 0, 
only the states m — ±S are appreciably populated and 
each molecular spin becomes effectively a two-level sys- 
tem ("spin-up" and "spin-down" states). In this "Ising" 
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FIG. 5: Temperature dependence of the equilibrium linear 
and nonlinear susceptibilities normalised by their values for 
isotropic spins (Brillouin limit). Upper panel: linear sus- 
ceptibility; lower panel: nonlinear susceptibility. Results are 
shown for a single crystal of Mni2 with the field parallel to 
the anisotropy axis o and a powder sample •. The lines are 
theoretical results for classical spins (dotted) and quantum 
spins (solid). 



limit, xt and X3T can also be calculated explicitly 

Xlsing 



( gM B) 2 S 2 2 
JVA ; — TZ, COS 1f> 



k B T 



X3Is 



(9) 
(10) 



Here we see that the bare 1/T and 1/T 3 dependences are 
recovered, so that the normalised susceptibilities become 
again temperature independent. Comparison of these ex- 
pressions with Eqs. JSJ-JHl reveals that, for a single crys- 
tal at tp = 0, xt/xiso and X3T/x3iso should increase, re- 
spectively, by an overall factor of 2.7 [= 3 x S 2 /S(S + 1)] 
and 12.3 (= 15 x S 4 /S(S + 1)[S(S + 1) + ±]), when de- 
creasing temperature. The experimental curves approach 
indeed these values at low temperatures (T < 10 K in 
Fig. [SJ although they become even larger. This is prob- 
ably due to the interactions which, as we have seen, en- 



hance the magnetic response at low temperatures. 

Anyhow, these results suggest the measurement of the 
temperature-dependent reduced susceptibilities as a suit- 
able tool to estimate the anisotropy parameters of su- 
perparamagnets. In this respect, the advantage of the 
nonlinear susceptibility becomes evident when dealing 
with systems with randomly oriented axes. Then the 
ratios between the "Ising" and isotropic limits decrease 
significantly [cf. Eqs. ©-© with Eqs. Pl-flltty]. Indeed, 
Xt/Xiso becomes nearly T-independent, whereas the re- 
duced nonlinear susceptibility still retains a sizable vari- 
ation with T (by a factor ~ 2.5). This is experimentally 
confirme d by measurements on a polycrystalline sample 
(Fig. |SJ| 33]. Thus we see that, contrary to the linear 
response, X3T keeps information on the anisotropy even 
for superparamagnets with axes distributed at random. 
This is the case most often encountered in nanoparticle 
systems [lj but also for single-molecule magnets when 
deposited on surfaces [34[ or inside porous materials [35| . 

The considerations above can be supported by direct 
diagonalization of the Hamiltonian J7J) . The results (solid 
lines in Fig. |SJ) exhibit the same trends as the experi- 
ments, both for parallel axes and after averaging over 
random orientations. Full agreement is precluded by the 
effect of interactions, on the low-T side, and by the pop- 
ulation of excited multiplcts with S ^ 10 at high T, as 
discussed above. 

Before concluding this section, there is an additional 
feature that deserves to be commented upon. Consider 
the theoretical behavior of X3T in classical spins, also 
plotted in Fig. [3(|. We see that, although classical 
and quantum calculations predict the crossover from the 
isotropic to the Ising limits, the classical susceptibilities 
are shifted towards lower temperatures. This shift can 
be seen as a manifestation of the quantum, discrete na- 
ture of the energy spectrum of Mni2 . The finite energy 
gap between the two lowest quantum levels, f^o, leads to 
a faster (exponential in D/T) convergence to the Ising 
limit, whereas classically this limit is only approached 
with a slow power law in D/T (see Appendix I A 2fl . 



IV. 



DYNAMICAL SUSCEPTIBILITIES 



In this section we turn our attention from the equi- 
librium to the dynamical response. We begin reviewing 
briefly the behavior of the nonlinear susceptibility X3 in 
the classical case. This allows to introduce some basic ex- 
pressions valid also for quantum superparamagnets; then 
we present the experimental results for Mni2. 



A. Classical superparamagnets and modelization 

The dynamical nonlinear susceptibility of classical 
spins was theoretically found to be very large and, in con- 
trast to the linear susceptibility, non-trivially sensitive to 
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the spin-bath coupling strength A 0, |3tJ . The "damp- 
ing" A measures the relative importance of the relaxation 
and Hamiltonian (precessional) terms in the dynamical 
equations 0, . Thus 1/ A is of the order of the num- 
ber of precessional turns that the spin executes in the 
deterministic spiraling down to the energy minima. 

The contributions to the nonlinear response of the lon- 
gitudinal and transverse components of the field are cap- 
tured by a simple formula involving the low-H expansion 
coefficients and gj_ of the relaxation rate [3^ 



r~r [i + i( fl ||£f + g±e±) 



(ii) 



where Tq = T\h=q and £ = g^SH /k^T The ex- 

pression for the nonlinear susceptibility oscillating with 
the third harmonic of the field |37| can be found in Ap- 
pendix It is easy to find the counterpart for the \3 
oscillating with e as used in this work, which reads 



(ffMB) 4 ^ 4 
IUJT 



X3 = 



= -N t 



COS Ip 



IUJT 



2(1 



IUJT) 



(<7|| cos 4 tp + g± cos 2 ip sin 2 ip) 



(12) 



Here the Ising approximation for the equilibrium parts 
has been used (this works fine at temperatures around 
the blocking temperature ~ 5K; see Sec. 11111 . The 
longitudinal part, proportional to cos 4 "ip, is maximum at 
tp = (in absolute value); the "transverse" contribution 
associated to g± becomes zero both at tp = and 7r/2, be- 
ing maximum at ip = it /A. Equation Ijl2|) shows that the 
magnitude, signs, and the u> and angular dependences of 
X3 are determined by the competition between the rate 
expansion coefficients gu and g±. Therefore, measure- 
ments of those dependences can provide valuable infor- 
mation on the different contributions to the spin reversal. 

For classical superparamagnets, where only the ther- 
moactivation operates, the rate expansion coefficients are 
given in the considered low-temperature range by |37ll38| 



= 1 , 9X = F(X)/2 , 



(13) 



where F > is a function of A (and T) . For strong damp- 
ing, F — ► 1, so that gu and g± are of the same order. In 
contrast, in the weak damping regime (governed by the 
precession), one has F oc 1/A and g± becomes very large, 
dominating the nonlinear response. Then one can find 
that the real part X3 ^> \x3t\ but with x'z/Xzt < 
(i.e., the sign is reversed with respect to the equilib- 
rium value) . This phenomenology is equivalent to that of 
the third harmonic nonlinear susceptibility |l7l l3^ | , since 
both quantities exhibit similar dependences on A and to. 

It is worth mentioning that in the derivation of Eq. I|12|) 
rather general assumptions are invoked (Appendix I A 3() . 
Therefore, the functional form obtained is quite generic 
and valid for classical as well as quantum superparam- 
agnets. Then, the relevant information on the quantum 
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FIG. 6: Nonlinear susceptibility vs. frequency at T — 5 K. 
Results are shown for different angles tp between the applied 
field and the anisotropy axis. The data are normalized by 
the equilibrium X3T measured at ip = 0. The relaxation time 
r was obtained from ^-dependent linear susceptibility mea- 
surements (as those of Fig- El - Upper panel: real part X3i 
lower panel: imaginary part x'i- The lines are obtained using 
Eq. 1121 with gu — —260 and g± — 0. The inset shows the 
sharply different classical prediction for %p = 0. 

reversal mechanisms will be incorporated by the rate ex- 
pansion coefficients gu and g±. Naturally, they could be 
very different from their classical counterparts 1131) . 



B. Nonlinear dynamical susceptibility of Mni2 

After these theoretical considerations let us turn to the 
experiments on quantum nanomagnets. Figure[S]displays 
frequency-dependent measurements of the nonlinear sus- 
ceptibility of Mni2 performed at constant T for differ- 
ent angles tp. They demonstrate that the result already 
shown in Fig. namely, that x'3 becomes much larger 
than its equilibrium value \3T near the blocking tem- 
perature, is a dynamical effect not caused by magnetic 
ordering or some kind of "freezing" (interactions also en- 
hance the susceptibilities, but by a much smaller factor). 
We know that classically one can also have |% 3 | ^> |x3t|, 
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but here X3/X3T > 0, that is, the peak of the measured 
nonlinear dynamic susceptibility is reversed with respect 
to the classical prediction. 

From Eq. I|12|) we see that at tp = there is no con- 
tribution of <7j_ to X3- I n addition, the first term in that 
equation, which has a Debye-type profile, cannot provide 
|X3 1 > IX3t| because Re[x/(1 + ix)] < x- Therefore, 
the maximum observed in Fig. should be due to the 
<7|i contribution. There is a result relating the height 
of the susceptibility peak Xalmax with the combination 
Q(ip) = cos 2 ip + g± sin 2 tp of the relaxation rate coef- 
ficients ( Appendix IA 3|) 

XsUx/XstI^ - -cos^ W)/4 . (14) 

Therefore, the positive sign of the maximum of 
X 3 (u;)/x3T at ip = entails Q < 0. But Q\^=o = g\\, en- 
tailing that the relaxation time t = 1/T becomes longer 
as increases. No classical mechanism can account for 
this; actually, g\\ = 1 in the classical model, which gives a 
Xs/X3t\4>=o smaller than 1 and decreasing with u> (inset 
of Fig. 0), in sharp contrast to the measured xs- On the 
other hand, it is well-established [21II2II2I I2I I2I 
that in Mni 2 the suppression of tunneling by a longi- 
tudinal field strongly reduces the relaxation rate (as it 
breaks the degeneracy between initial "spin-up" and fi- 
nal "spin-down" stateSjinhibiting the tunnel channels). 
As suggested in Ref. this effect provides the g\\ re- 
quired, both negative and large, to account for the ex- 
perimental X3 °f Mn 12 . Thus, we see that the known 
field-suppression of tunneling shows up in the nonlinear 
response as a distinctly quantum contribution, with its 
sign reversed with respect to the classical case. 

When ip > 0, the detuning quantum contribution co- 
exists with the transverse g± contribution. Still, the data 
of Fig. suggest that X3 cos4 ip holds approximately. 
We can check this by accounting again for Eq. H14fl and 
plotting the maximum of x'3 vs - cos 4 ip (Fig. [7J. This 
yields an almost straight line indicating that in Mni2 the 
coefficient g\\ overwhelmingly dominates gj_, which in the 
classical case embodied the precessional contribution and 
could be sizable (we return to this issue in Sec. IVjl . 

Thus X3(w,^) provides direct experimental access to 
the relaxation-rate expansion coefficients g\\ and g±, 
which contain information on the spin relaxation mecha- 
nisms. Besides, from the sign of the X3( w ) peaks we can 
infer whether the spin reversal is dominated by a classical 
mechanism or by quantum processes. The consistency of 
our analysis can be ascertained by comparing the exper- 
imentally determined rate T (obtained from Debye fits 
of x)> with the rate reconstructed from the expansion 
T = r (l + Q£ 2 /2 + ...), using the Q\^ = o extracted 
from the x'3 maxima via Eq. (|14|) (QMn 12 \i>=o — —260). 
Figure |H1 shows the good agreement between both results 
in the weak field regime, supporting our interpretation. 
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FIG. 7: Evolution of the maximum of x'ai 10 ) with the angle 
ip of the applied field at T = 5K. The symbols are the exper- 
imental results. The dashed lines are obtained from Eq. (I12B 
using the experimentally determined — —260 and the clas- 
sical transverse contribution g± = F(A)/2 computed for dif- 
ferent values of the phenomenological damping constant A. 
The solid line corresponds to theoretical calculations with a 
Pauli quantum master equation (see the text). 
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FIG. 8: Relaxation rate of Mni2 at T = 5K as a function of 
the field applied along the anisotropy axis (ip = 0). The sym- 
bols are the T's obtained from Debye fits of x( w > H). The lines 
are calculated as V = ro(l+Q£ 2 /2), where £ = gfis SH^ /kBT 
and Q comes from the maximum of x'si 10 ) v i a Eq. 11411 . The 
classical prediction (Q = 1) is shown for comparison. 

V. DISSIPATION VS. DECOHERENCE IN 
SINGLE-MOLECULE MAGNETS 

We have seen in the previous section that the trans- 
verse contribution to the nonlinear susceptibility is nearly 
absent in Mni2 (indeed the experiments are consistent 
with g± = 0). Classically, g± incorporates the preces- 
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sional contribution, which can be large for weak damp- 
ing a<i mill. Still, as its long t indicates, Mn 12 is 

expected to be quite "underdamped" (in the sense of en- 
ergy dissipation). Thus, it seems that some other process 
makes the spin to loss its intrinsic, or coherent, preces- 
sional dynamics and appear, when seen through X3(w), as 
"overdamped" . In this last section we will try to recon- 
cile these results (gj_ ~ and long To) invoking an effect 
of the coupling to the bath absent in classical physics — 
decoherence. 



A. Experimental bound for the effective damping 

Before starting let us quantify a lower bound for an ef- 
fective A of Mni 2 . To this end we generate X3( w ) curves 
using Eq. i|12|) with the <?|| experimentally determined 
from X3( w )|i/j=0i while assuming g± = F{\)/2 as in clas- 
sical superparamagnets. In this way we compute x^lmax 
vs. cos 4 ip for several A and compare with the experimen- 
tal results (Fig. CJ). The best agreement is obtained for 
large A, which in fact yields small g± and hence almost no 
cos 2 tp sin 2 tp contribution. However, taking into account 
experimental uncertainties as well as the smaller sensi- 
tivity of g± to large A, gives a lower bound of A > 0.01. 

In the classical equations of motion the damping pa- 
rameter A is a measure of the relative weights of their 
relaxation and precessional terms. It can then be ex- 
pressed as the Larmor precession period tl = 2n/"fH aj 
(in the anisotropy field H a ; 7 = g^e/ft) divided by some 
timescale of relaxation. For the latter we can use the 
classical prefactor in the Arrhenius law 0, Ei2 



To = 



(15) 



where the reduced anisotropy barrier is a = U^/k^T. In 
our experiments at T = 5 K we have a ~ 14, whence 



A ~ 0.04 x tl/to . 



(16) 



Experimentally To ~ 3 x 10~ 8 s in Mni2. Its anisotropy 
field can be obtained from magnetization measurements 
(along the hard plane) or from the ground-state transi- 
tion frequency Qo/h [Eq. ©], getting tl ~ 3-4 x 10 -12 s. 
This gives tl/t ~ 10~ 4 and hence an effective A many 
orders of magnitude smaller than the lower bound A > 
0.01 extracted from X3- 

The estimation of tl/t is in agreement with the men- 
tioned underdamped character of Mni2 (in the sense of 
long-lived levels). However, in case of g± having some 
precession- type contribution akin to F(X), the above 
A oc tl/to would not be the parameter entering there. 
Otherwise a giant gj_ tx 1/A would dominate X3 (leading 
to X3/X3T < and Xalmax not proportional to cos 4 tp), 
which is clearly not seen in the experiments. Still, 
it might be that, instead of the damping, some other 
timescale limits the precession; as this is a type of coher- 
ent dynamics, we can call such time a decoherence time 



Tdoc- This should replace t in the effective A cx t^/t^cc- 
To get a A consistent with the lower bound A > 0.01 ob- 
tained from X3i the time Tdec should be far shorter than 
t ~ 10~ 8 s, actually T dcc ~ 1-1.5 x 10~ n s. 



B. Approximate quantum treatment 

A critical assessment of the considerations above is 
in order before proceeding. They have been based on 
stretching the classical idea that g± should include some 
precessional-type contribution. Further, they assume 
that this contribution is controlled by a parameter re- 
lating the Larmor period with some scale limiting the 
time allowed to the spin to precess (either by damping 
or by some loss of coherence). By definition, however, 
g± accounts for the effects of H± on the relaxation rate 
[Eq. JTIJ]; besides, X3( w ) gives direct access to g\\ and 
g±. Then, one would expect that a pure quantum ap- 
proach for r, including the bath coupling and coherent 
dynamics (tunnel and precession) could account for the 
experimental X3 without recourse to classically precon- 
ceived notions. 

An exact quantum treatment, unfortunately, is diffi- 
cult because one must deal with the full density-matrix 
equation including the intrinsic (Hamiltonian) dynamics 
plus the effects of the bath (damping and decoherence). 
However, as this can be handled in various limit cases, 
we shall attempt a discussion based on the correspond- 
ing partial solutions for the quantum T, with the hope of 
shedding some light on the physical origin of the results. 

The dominant terms that enter the relaxation rate de- 
scribing quantum tunneling via a pair of nearly degener- 
ate states |±m) has a Lorentzian shape as a function of 
the longitudinal bias £ m = 2mgii^H^ [40| 



2r r 



'£2 1 a/ +h 2 T 2 



exp(-C/ m /fc B T) . (17) 



Here T m is the probability of decaying to other levels via 
the absorption or emission of phonons (i.e., it is approxi- 
mately 1/to), A m is the tunnel splitting of the pair |±m) 
induced by terms not commuting with S z in the Hamil- 
tonian, and U m is the energy of the levels. The width 
of the Lorentzian introduced |4l|, w 2 m = A 2 m + tPT^, 
interpolates between the results that can be obtained in 
the two limit cases of (i) large coupling H/tq 3> A m [24|. 
where w m ~ h/ro, and (ii) weak coupling H/tq <C A m 
[25|. in which w m ~ A TO . 

Performing the second _ffj^-derivative of Eq. I|17|l gives 
the corresponding g±. For the parameters of Mni2 it 
results that the main contribution at zero longitudinal 
bias comes from the derivative of the quotient A^/w^ 



1 9 2 r 



9-L 



T dH 2 



d 2 



(A m /w m ) dB\ 



{A m /w m ) . (18) 
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Now, for large coupling h/r ^> A m , we have w m — h/r . 
Then the relaxation rate depends on the ratio between 
To and the tunneling time H/A m . A transverse field eases 
the tunneling and significantly increases A m , and hence 
r. As r is then quite sensitive to H±, we can have large 
gx, in analogy with the classical situation. On the con- 
trary, when H/tq -C A m , we have w m ~ A m and hence 
^m/ w m — 1- Then the rate becomes quite insensitive to 
A m (and hence to H±), leading to a small g±. 

In Mni2, where H/tq ~ 0.2 mK, both situations are 
possible. The reason is the exponential increase of the 
tunnel splitting A m with decreasing |m|, going from the 
sub-nanoKelvin regime for the ground levels m = ±S 
to some tenths of K for \m\ < 2. Therefore, the rela- 
tion between H/tq and A m depends on which tunneling 
path (i.e., which pair ±m) gives the dominant contri- 
bution to r. If tunneling proceeds via the deep levels, 
where H/tq S> A m , we would find a large g±. In con- 
trast, when tunneling occurs through the excited levels, 
one has H/tq <C A m and hence small gj_. 

At this point it is important to bring into the discussion 
the effect of environmental bias fields (due to intermolec- 
ular dipolar interactions or the hyperfine interaction with 
the nuclear spins of the Mn ions). They produce a distri- 
bution of bias £ whose typical width is of a few tenths of 
K (of the order of the Curie- Weiss 9). These bias fields 
enter as A^/(£^ + w^) in the rate expression l(T7)l. re- 
placing the bare A m /w m and supressing tunneling when 
w m <C Cm (Fig.[TJ. Taking into account the order of mag- 
nitude of £ m , the bias effectively blocks tunneling via the 
large \m\ (deep) channels, those that would provide large 
gx- Tunneling becomes possible only for A m > £ m , but 
for those upper levels H/tq <C A m , leading to small g±, 
in agreement with the experiments. 

We can support this picture with direct numerical 
calculations. An approximate (Pauli) quantum master 
equation, which works well when tunneling occurs un- 
der weak damping conditions and that incorporates the 
effects of environmental bias fields, was used to study sev- 
eral problems in Mni2 pH l42l |43| . We have implemented 
it to address the nonlinear response problem, mimicking 
the experimental protocol and calculating \' an d x" vs - 
H. Results are shown in Fig. (solid line; see also Fig. 3 
of Ref. They account well for the measured non- 

linear susceptibility; in particular, for the nearly cos ip 
dependence of X3 associated to small g±. 

We would like to also provide a physical picture in the 
limit cases discussed (H/tq much larger or smaller than 
A m ). To this end, let us discuss the total energy of the 
spin plus the bath, treating their interaction perturba- 
tively. Then, time-dependent perturbation theory leads 
to the celebrated time-energy "uncertainty" relation. In 
particular, for t > h/ AE the dominant transitions are 
those conserving the total energy. On the other hand, for 
a spin in a |m) state, which is not an exact eigenstate of 
the Hamiltonian, the energy uncertainty is of the order of 
the tunnel splitting A m . Well, consider now a transition 
to such an |m) state; although the spin could at short 



times remain there, for times longer than r$ = H/A m 
it will have to become an energy eigenstate. Then the 
wavefunction consists of a superposition of spin-up and 
spin-down states | ±m) , delocalized between both sides of 
the barrier. 

We can now revisit the limit cases discussed above. 
Consider first the strong-coupling case H/tq 3> A m . For 
times shorter than the decay To, the time-energy un- 
certainty allows the existence of superpositions of en- 
ergy eigenstates, which can be localized on either side 
of the barrier (~ |±m)). These wavepackets may exhibit 
Hamiltonian dynamics including tunnel oscillations and 
precession. Under these conditions the rate T is quite 
sensitive to A TO , as it controls the probability for the 
spin to have tunneled before a time To- It is this sensi- 
tivity to A m , and in turn to H± , which can lead to large 

g±- 

When, by constrast, H/tq <C A m the "semiclassical" 
wavepacket around \m) delocalizes in the tunneling time 
t$ = h/ A m , evolving towards an energy eigenstate due 
to the uncertainty principle. Then, there is no coherent 
oscillation between |m) and \—m), neither precession of 
the (averaged) transverse components, since this is a sta- 
tionary state. The dependence on A m (and hence on H±) 
is then minimized, as the wavefunction is already delocal- 
ized between the spin-up and spin-down states, leading to 
small g± values. Note that under these conditions the co- 
herent (precessional) dynamics is not limited by the level 
lifetime To but by the much shorter "decoherence" time 
to attain a diagonal density matrix. Therefore, the Td ec 
introduced heuristically above can be identified with this 
t$ = h/A m . For \m\ = 2, 4, we have A m ~ 0.7-0.02 K, 
which give t$ ~ 10~ n -4 x 10~ 10 s. These values are con- 
sistent with the estimation, based on Xai^), of the Td ec 
required to get A ~ 0.01 (which yielded Td oc ~ 10~ n s) 



VI. SUMMARY AND CONCLUSIONS 

The single-molecule magnet Mni2 is a model system 
for the study of thermal-equilibrium properties of spins 
with magnetic anisotropy, as well as the dynamics of a 
quantum mesoscopic system subjected to the effects of 
a dissipative environment. In this system, to the rich 
physics of classical supeparamagnets, quantum effects are 
incorporated. In this article we have investigated exper- 
imentally the equilibrium and dynamical nonlinear re- 
sponses of Mni2- The nonlinear susceptibility X3 was 
underexploited in this field in spite of having, when com- 
pared to the linear susceptibility, an enhanced sensitivity 
to several important characteristics of the system. 

We have shown the sensitivity of the equilibrium X3 to 
the magnetic anisotropy parameters of the spin Hamil- 
tonian. As in the classical case, the anisotropy leads 
to an extra temperature dependence of X3 that, in con- 
trast to the linear susceptibility, persists for randomly 
distributed anisotropy axes. Therefore, the measurement 
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of X3 (T) can be exploited to estimate the anisotropy con- 
stants even in powdered samples or in systems deposited 
on surfaces or in porous materials. 

The analysis of the dynamical \3, with help from a 
generic but simple formula [Eq. provides valu- 

able information on the intrinsic dynamics of the system. 
Specifically, X3( W ;V0 gives access to the relaxation-rate 
field-expansion coefficients g\\ and g± [Eq. 1)11(1 ]. which 
contain information on the mechanisms of spin reversal. 
The experimental nonlinear response of Mni2 is found to 
be consistent with the established scenario of thermally 
activated tunnel via excited states. Thus, the strong de- 
crease of the relaxation rate due to the (longitudinal) 
field-detuning of tunnel levels manifests itself in \3 as a 
distinct quantum contribution, with a sign opposite to 
the classical case. Then, from the signs of the %3 vs. u> 
peaks one can estimate if quantum effects play a role in 
the dynamics of the studied nanomagnet. Finally, from 
the analysis of the angular dependence of X3(w) we have 
estimated a bound for the decoherence time required to 
attain a diagonal density matrix due to the phonon-bath 
coupling. The so obtained r$ is much shorter than the 
level lifetime To and is the responsible for a fast loss of 
coherent dynamics like tunnel oscillations or precession. 
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APPENDIX A: ANALYTICAL EXPRESSIONS 
FOR VARIOUS QUANTITIES 

In this appendix we derive a number of analytical for- 
mulas used in the discussions of the main text. First, 
the equilibrium nonlinear susceptibility of a spin 5 in 
the isotropic limit (from the Brillouin magnetization). 
Then, corrections due to finite magnetic anisotropy to 
the equilibrium linear and nonlinear susceptibilities in 
the opposite Ising limit. Finally, we derive the frequency 
dependent nonlinear susceptibility 1)12(1 : we also analyze 
the zeroes, extrema, and signs of the formula for Xii 1 ^)- 
For simplicity, we omit throughout the appendix unessen- 
tial constants like Na, 5Mb, ^b, etc. 



1. Equilibrium isotropic susceptibilities from the 
Brillouin magnetization 

For isotropic spins, we can get the equilibrium linear 
and nonlinear susceptibilities by expanding the Brillouin 
function around zero field. The isotropic magnetization 



can be written as (y = H/T) 

M z = a cth (ay) — b cth (by) , a — S + ^ , b = i . 

Then, from the first terms of the small x expansion of 
the hyperbolic cotangent, namely cth (x) ~ 1/x + x/3 — 
x 3 /45, we get 

M ~ W„2 h 2\„, 1 i„A U A\„,3 4£ f „ ej , „ _u3 



'z — 3 



K - b 2 ) y - A (a 4 - & 4 ) y A = Xt H + X3T H 



Using nowa 2 -6 2 = 5(5+1) and a 2 + b 2 = 5(5+l) + i, 
we finally obtain 

5(5 + 1) 5(5+l)[5(5 + l) + ±] 

XT = ™ , X3T — 



3T 



45 T 3 



(Al) 

The first is the celebrated Curie Law and the latter its 
sought generalization for the nonlinear response. Note 
that there is not second order \2 since M z is odd in H . 



2. Finite-anisotropy corrections to the equilibrium 
Ising susceptibilities 

Now we go into the corrections due to finite anisotropy 
to the susceptibilities in the opposite Ising limit. In this 
regime, for uniaxial anisotropy, only the states m = ±5 
are populated at zero field, and each cluster becomes 
effectively a two-level system. We will consider the ef- 
fects of finite population of the first excited levels m — 
±(5 — 1). For simplicity, and to compare with known 
results in the classical case we assume the simplest 
form for the magnetic anisotropy H = —DS\ [cf. Eq. Q]. 

We start from the statistical-mechanical expressions 
for the susceptibilities at zero field [T^ |: 



XT = 



T 



X3T 



(si) -s (s* z y 

6T 3 



(A2) 



To get the equilibrium averages (S z ), we need first 
the partition function. Including the contributions of 
the ground states, m = ±5, and first excited states, 
m — i(5 — 1), we simply have 



Z=± e-/- 

m— — S 



2(e- £s / T + e- £s - 1 / T ) , (A3) 



where we have taken into account the zero-field degener- 
acy e- m = e m (yielding the factor 2). For e TO = —Dm 2 
the separation between adjacent levels is £ m — e m _i = 
— D(2m— 1). Thus, for the ground-state splitting one has 
£s — £s-i — —D(2S — 1) = — Oo, which corresponds to 
Eq. JSJ. Then, extracting a factor e~ £s ^ T in Z, we find 
the approximate partition function for our problem 



2e" CT (l 



-Qo/T) 



(A4) 



where we have also introduced the reduced anisotropy 
barrier a = DS 2 /T. 
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Next, we compute the moments (S k ) oc ~J2 rn m k e e ™/ T 
required in Eqs. (|A2() within the same approximation: 



(s k z ) 



1 



S- l\ k 



s 



-fio/T 



(A5) 



We have considered even k, yielding the factor 2; oth- 
erwise the moments vanish (as they are computed at 
H = 0). Then, dividing by Z as given by Eq. I|A4J) . 
we have 



(S k M ) = S> 



1 



(A6) 



l + e -n /T 

Inserting these moments for k = 2, 4 in Eqs. (|A2|) . we get 
S 2 1 



XT 



T 



1 + c- f VT 

+ W) 4 e 
6T 3 \ 1 + e- n °/T 



g4 p+^V^o/T 
X3T — < : tt-^ A 



fio/T 



l + e-^o/T 



which are the desired susceptibilities incorporating the 
effects of finite population of the first excited levels. 

The above results show that the corrections to the Ising 
limits xt = S 2 /T and X3T = -S 4 /3T 3 are functionally 
exponential. This can be seen more explicitly as follows. 
Note first that (3 = e~°°/ T is a small parameter in the 
considered approximation (for Qq ~ 15 K and T ~ 5K, 
we have (3 ~ 5-1CT 2 ). Then, introducing f k = (1-1/S) k , 
the moments can be approximated binomially, (l + x) a = 
l+ax---, giving (S k ) = (l+/3/ fe )/(l+/3) ~ l-/3(l-/ fe ). 
To illustrate 1 - / 2 = (25 - l)/5 2 gives explicitly 



S 2 

XT = y 



1 - 



25-1 
^ 2_ 



-n /T 



(A7) 



and a result structurally similar for X3T- These exponen- 
tial corrections are to be compared with the pow er law 
corrections in the classical asymptotic results 



XT 



El 

T 



1 



X3T 



s*_ 

3T 3 



(-§) 



Thus, while in the classical case the corrections enter as 
inverse powers of a = DS 2 /T, for quantum spins they are 
exponential in D/T, leading to a much faster approach 
into the Ising regime as the temperature is lowered. 



3. Dynamical nonlinear susceptibilities 

Finally, we proceed to derive Eq. I|12|) for the frequency 
dependent nonlinear susceptibility x^ used in this work. 
The ordinary nonlinear susceptibility, denoted here X3 j 
is defined in terms of the third harmonic of the response 
to an alternating field. For classical spins with uniaxial 
anisotropy a formula for x^t^ was derived |37l ] from a 



system of low T balance equations 01 obtained from 
Brown Fokker-Planck equation , which reads 



(3o 

X 3 



5 4 
3T 3 



3iwr (g\\b' 



l + 3iwr 2(1 + iwr)(l + 3iwr) 

(A8) 

Here 611 = H\\/H = cosip and b± = H±/H — sin?/; are 
the direction cosines of the probing field, while g\\ and 
g± are the expansion coefficients of the relaxation rate as 
given by Eq. (JTTJ), T ~ r [l + \{g\\b\ + <?±&iK 2 ], with 
To = r|w =0 and £ = SH/T (these g's differ from those of 
Ref. [33 by a factor 1/2). Accepting the validity of the 
balance equations for the number of spins pointing up or 
down, the derivation of Ref. J37J leading to Eq. QASj l is 
also valid for quantum spins (with uniaxial anisotropy). 

An alternative nonlinear susceptibility [ToL ITU Il4| , the 
one we use in the main text, can be obtained from the first 
harmonic of the response in the presence of a weak static 
Afield if as X (u,H) = X i{u) + Zxf ] H 2 + . . . [cf. Eq. Q]. 



jThen, Xz^ can be obtained from the first harmonic as 



d%x\ 



3! x 3 



(A9) 



Experimentally, as discussed in Sec. [HJ one can mea- 
sure x f° r several H and get x^ by polynomial fitting. 
Both x^^ an d X3 coincide in the limit u — > with 
the thermal equilibrium X3T- In addition, as we will see 
here, they have qualitatively similar dependences on the 
damping, ui, T, angle, etc. However, Xa^ presents the 
experimental advantage of not requiring high-harmonic 
detection and processing. 

To derive a formula for ^3 one can proceed as in 
Ref. using a system of low T balance equations. Here, 
however, we shall present a more direct derivation start- 
ing from the Debye form for the first harmonic re- 
sponse. Let us write this as xi^i H) — xs + Ax/1 + i lot, 
where Ax = Xt — Xs- This one-mode relaxation form 
is valid in weak enough fields just account- 

ing for the H dependences of the x' s an d t. We can 
take advantage of this to get Xa by simple differenti- 
ation x^ — &hX\o/3" In addition, in the low tem- 
perature regime we approximate the susceptibilities by 
their Ising limits (this is supported by the experiments 
of Sec. IIIII . Consistently, the magnetization is given by 
M z /S = 6||th(£||), with ^ = b\\£. Then, we write the 
Debye law as 



1 



A x d H M z 



(A10) 



This low T approximation entails to disregard xs ^ Xt~ 
Xs and d'jjXs *C X3T H2, and corresponds to the 2-state 
approximation in the balance equations approach |37j . 
Based on the above considerations we simply write 



X = d H M z D, D 



(All) 
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and proceed to differentiate to get Xs = c^xlo/3!, using 
eventually T = T Q (1 + Qt; 2 /2). To work out the second 
derivative we use the "binomial" formula (/ g)" = f"g + 
2f'g' + f g" , which results in 

d 2 H x = d 3 H M z ■ D + 2 d 2 H M z d H D + d H M z d 2 H D . (A12) 

Now, using the evenness of the rate on H, we have 
d?/r|o — 0, which has two important consequences 



equilibrium value X3 — Xa /X3T' 



X3 



l + ix 2 (1 + ia;) 2 



i = wt, (A15) 



where Q = Q/b 2 = g\\ + g±(b±/b^) 2 . Multiplying now by 
the conjugate denominators we readily separate the real 
and imaginary parts 



d H D =0 



Wo = 



iu;d 2 H T 



(r- 



1W 



Then, from T = r (l + Q£ 2 /2) and d H = (S/T)% along 
with 3hM z |o = Xt and djjM z \o — 3!x3T, we finally ob- 
tain 



X3 



X3T 
1 + \LOT 



+ 



XtS 2 



\LOT Q 



T 2 6(1 + icjr)(l + iwr) 



(A13) 



where now r = 1/ro is the zero-field relaxation time. 

To compare with Eq. I|A8(1 for xf^i we recall that 
M z /S = 6||th(£||) and use the expansion th£ — £ — ^£, 3 , 
whence d H M z = (S 2 /T)(bj - b^ 2 ) and d 3 H M z \ = 
6(-S" 4 6|/3T 3 ). Then, introducing the explicit expression 
for Q = g\\b 2 + g±b 2 ± , we arrive at [cf. Eq. JTSJ)] 



X3 



~3T 3 



iujT (g\\bf. 



g±b 2 b 2 ± r 



1 + iujT 2(1 + iwr)(l + iut) 

> A 1 I ) 

This result is generic and valid for classical and quan- 
tum spins, inasmuch as the starting Debye x provides an 
adequate description. In general, g\\ and g± will not be 
given by the classical result JT3J, but they incorporate 
quantum contributions to the relaxation rate. 

The expression derived for X3 shows a close struc- 
tural analogy with that for the third harmonic nonlinear 
susceptibility [cf. Eq. I|A14|) with (|A8|) ]. Indeed, we have 
kept the factor (1 + iuiT) 2 without squaring in Eq. (1 A 14(1 
to enhance the analogy; replacing lot — > 3wr both quanti- 
ties exhibit almost the same frequency dependece. Com- 
parison of Eq. I|A14|) with l|A8(l also reveals similar de- 
pendences on T, (?|| and g±, as well as on the angle ip. 
This supports our repeated claim about the analogous 
qualitative dependences of xi^ and ^3 and, in turn, 
our choice of the first harmonic response on the basis of 
its experimental convenience. 

We conclude this appendix finding extrema and zeroes 
of X3(w). This will help to explote having an analytical 
expression for the nonlinear susceptibility when analyz- 
ing experiments. First we normalize Eq. l|A14f) by the 



X 3 = 



1 + x 2 (l + x 2 ) 2 ' 



-Xs 



x Q x(l — x 2 ) 
l + x 2+1 2 (l+x 2 ) 2 



Let us first compute where the imaginary part crosses 
the lot axis (see Fig. bottom panel). To find this re- 



quires to solve x(l + x-) + (Q/2)x(l 
from x — 0, one finds the following zero: 



Q 



IQI»i 



1 + 



Q 



0. Appart 



(A16) 



To get the large Q approximation we have used the bi- 
nomial formula (1 + x) a — 1 + ax ■ ■ ■ twice (to work the 
denominator and then to take the square root). Well, 
note that for \Q\ < 2 one has x 2 < and hence the imag- 
inary part does not cross the x axis. However, when one 
finds the crossing, the following criterion holds: if x z > 1, 
then Q > 0, whereas for x z < 1, one has Q < 0. Thus, 
inspection of the crossing of the lot axis (below or above 
lot = 1), could provide information on the classical or 
quantum character of the spin-reversal dynamics. 

Let us now find the extrema x m of the real part x 3 - 
Writing x' 3 = [1 + (1 — Q)x 2 }/(1 + x 2 ) 2 one sees that 
dx' z /Ax = implies — (1 + Q) + (1 — Q) x 2 n — 0, whence 



Q + i 

g-i 



— > 



1 

Q 



(A17) 



There is a maximum, or minimum, when \Q\ > 1: in such 
case, x m > 1 entails Q > 0, whereas x m < 1 yields Q < 0. 
We conclude giving the value of x'z at its extremum 



XsOm) = - 



4Q 



IQI»i 



(A18) 



which is the result used in Sec. |IV| to get g\\ from the 
maximum of X3 at ip = (then Q = g^). Note finally that 
the signs of the x'3 peak and of Q are always opposite, 
which can help ascertaining the quantum contribution to 
the nonlinear dynamics. 
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